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Studies of free particles in low-dimensional quantum systems such as two-leg ladders provide 
insight into the influence of statistics on collective behaviour. The behaviours of bosons and fermions 
are well understood, but two-dimensional systems also admit excitations with alternative statistics 
known as anyons. Numerical analysis of hard-core Z 3 anyons on the ladder reveals qualitatively 
distinct behaviour, including a novel phase transition associated with crystallisation of hole degrees 
of freedom into a periodic foam. Qualitative predictions are extrapolated for all Abelian Z, anyon 
models. 


The study of low-dimensional quantum systems 
through both theory and experiment provides a unique 
window into the collective behaviours of quantum many- 
body systems [1-3]. Competition between interaction 
terms, hopping terms, filling fractions, and geometrical 
frustration may give rise to an immensely broad range 
of different phases of matter, and transitions may be ob¬ 
served between regimes which are insulating, conductive, 
or superconductive, ordered or disordered, and many 
more. Moreover, the statistics of the particles in a model 
may significantly affect its phase diagram [e.g. 4-9]. 

Other than the ID chain, one of the simplest low¬ 
dimensional quantum systems is that of a hard-core non¬ 
interacting gas on a two-leg ladder. An infinitely long 
ladder may be considered to be a quasi-lD system, al¬ 
though at incomplete fillings it is possible for particles 
to braid around one another and thus the phase diagram 
is affected by particle statistics. The behaviours of both 
fermions and hard-core bosons on the two-leg ladder have 
been extensively studied, and are known to generate rich 
phase diagrams dependent both on the ratio of the rung 
and leg hopping terms and on the value of the chemi¬ 
cal potential [9]. However, these two examples are not 
exhaustive as there are other choices for particle statis¬ 
tics in one- and two-dimensional systems. These anyonic 
statistics are of fundamental interest as exotic and as- 
yet poorly-understood phases of matter. Efforts to un¬ 
derstand their behaviour are also increasingly driven by 
their roles in phenomena such as the fractional quantum 
Hall effect (FQHE), [10-18], and as exotic excitations in 
superconducting heterostructures and nanowires [19-28]. 

In this paper 1 consider a simple example of anyonic 
statistics, the Z 3 anyon model, which is realised in the 
FQHE at 1 / = I [11] and in coupled wire constructions 
[29, 30]. Whereas a bosonic system is unchanged under 
pair exchange and a fermionic system picks up a fac¬ 
tor of —1 when exchanging particles which are nearest 
neighbours with respect to some linearly-ordered basis, a 
system of anyons may undergo a unitary transformation. 
For Z 3 anyons, this transformation is a phase and is a 
third root of unity. 

The study of anyonic systems is computationally chal¬ 
lenging, as their non-trivial exchange statistics give rise 


to a form of the sign problem [31, 32]. However, the class 
of simulations known as tensor network algorithms are 
known to be relatively resilient against this effect [33-37] . 
Recent work on the incorporation of both Abelian and 
non-Abelian symmetries into tensor network algorithms 
[38-42] admits a generalisation to arbitrary anyonic sys¬ 
tems [43-45], and this has led to the recent development 
of anyonic tensor network algorithms such as anyonic 
Time-Evolving Block Decimation (TEBD) [46] and any¬ 
onic Density Matrix Renormalisation Group (DMRG) 
[47, 48]. In the present paper, the anyonic version of 
the infinite DMRG algorithm [48-51] is used to study 
Z 3 anyons on the two-leg infinite ladder. Performance is 
enhanced using precomputation techniques [41]. 

The Z 3 anyon model has two particle species, also 
known as charges, which are customarily denoted 1 and 2 , 
and a vacuum sector, denoted 0. Fusion of particles is 
equivalent to summation of charges, represented by addi¬ 
tion modulo 3 , 0 x 6 ^ {a+b) mod 3, and exchange of two 
neighbouring charges a and 6 in a linearly-ordered basis 
attracts a factor of or where i?®*” = exp( ^^j°^ ). 

Whether this factor is conjugated depends on which par¬ 
ticle passes behind, and which in front of, the other, with 
respect to some oriented, linearly-ordered basis. It fol¬ 
lows from these definitions that species 2 is the antipar¬ 
ticle of species 1 , with the model being invariant under 
simultaneous interchange of charge labels 1 and 2 and 
time reversal (complex conjugation). In this paper we 
consider only a single species of hard-core Z 3 particle on 
the ladder [52], with the resulting phase diagram being 
independent of our choice. For definiteness we may con¬ 
sider this to be species 1 . 

As with fermions we begin by imposing a linear order¬ 
ing on the sites of the ID ladder: 

(') 

c d 

This order is supplemented by an orientation—we spec¬ 
ify that this basis is viewed from the lower side of the 
diagram. Thus if a particle hops from site a to site b 
in the above diagram, it passes behind site c, and if one 
hops from c to d, it passes in front of site b. If a particle 
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FIG. 1. (i) Phase diagram for fermions on the 2-leg ladder, 
(ii) Phase diagram for hard-core bosons on the 2-leg ladder. 
Greyscale denotes filling fraction, from 0 (black) to 1 (white). 
Labels LLl and LL2 indicate 1-band and 2-band Luttinger 
liquids respectively in the fermionic system, SF denotes a su¬ 
perfluid phase for the hard-core bosons, and p = ... signifies 
an insulating phase with the specified filling fraction. The 
p = I phase is a band insulator for fermions, and a rung- 
Mott insulator for hard-core bosons. These results reproduce 
phase diagrams previously presented in Ref. 9. 




is also present on the site being hopped past, then this 
attracts a phase factor or . Using a graphical 

notation similar to that in Refs. 53 and 54 and the or¬ 
dering of Eq. (1), the Hamiltonian for hard-core particles 
with nearest neighbour hopping and a chemical potential 
may be written 



where a solid line indicates the presence of a particle 
(charge 1 ), the absence of a line indicates the vacuum 
sector (charge 0 ), and a grey line may be either. 

In simulating this system the accessory site variation 
of the iDMRG algorithm discussed in Ref. 48 provides 
exceptionally fast convergence towards the ground state, 
and comparison with Ref. 9 shows that for a refinement 
parameter D = 200, an accurate picture of the phase 
diagram for fermions and hard-core bosons emerges af¬ 
ter as few as 25 iterations, as shown in Fig. I, with the 
insulating and the Luttinger liquid or superfluid phases 
of these models being clearly visible. Application of the 
same technique to a system of Z 3 anyons yields the phase 
diagram shown in Fig. 2(i). Individual sections through 
this diagram were also computed with up to 300 itera¬ 
tions and D = 400 and were confirmed to show no quali¬ 
tative changes, with only minor quantitative changes on 
the order of one percent filling. The three-rung von Neu¬ 
mann entanglement entropy [Fig. 2(iv)] was com¬ 
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FIG. 2. (Golour online) (i) Phase diagram for hard-core Z 3 
anyons on the 2-leg ladder. Labels SFl and SF2 indicate su¬ 
perfluid regimes analogous to fermionic 1 -band and 2 -band 
Luttinger liquids, and p — . signifies an insulating phase 

with the specified filling fraction. HG labels the novel “hole 
crystal” phase described in the text, with thin black lines indi¬ 
cating the approximate boundaries of this phase, (ii) Section 
through diagram (i) at t±/t = 1 , showing (n) as a function 
of p/t. Results for Z 3 anyons are indicated by o, with results 
for fermions and hard-core bosons being extracted from Fig. 1 
and denoted -I- and x respectively, (iii) Section through dia¬ 
gram (i) at tx/t = 3. Key as for diagram (ii). (iv) Three-rung 
von Neumann entanglement entropy s[^'^ as a function of p/t 
for t±/t = 1. A discontinuity in at p/t ~ 2.04 shows the 
SFl/HC phase transition to be first-order. 


puted from the 3-rung reduced density matrix using 
the anyonic entanglement entropy formalisms of Refs. 55 
and 56. As Z 3 anyons have a quantum dimension of 1, 
these reduce to the usual expression 

S'f) =-Tr[,3(3)logp(3)]. (3) 

The most obvious feature of Fig. 2(i) is the loss of par¬ 
ticle/hole duality, corresponding to reflection symmetry 
of the phase diagram about p/t = 0. While this prop¬ 
erty holds for both fermions and hard-core bosons, its 
absence is unsurprising for Z 3 anyons: For both fermions 
and bosons this property may be identified with the ex¬ 
istence of a Z 2 symmetry associated with conservation of 
parity. Starting at the p/t = 0 half-filling state, any small 
perturbation introducing or removing a particle from the 
half-filled state causes interchange of the charge 0 and I 
sectors of the MRS chain, with a corresponding shift in 
the chemical potential being required to stabilise this 
change. For Z 3 anyons, this symmetry is broken. The 
dual of charge I is now charge 2, and addition or deletion 
of a particle on a site causes the charge sectors to incre- 
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merit or decrement by one, modulo 3. Thus the defects 
introduced in the MPS chain by insertion and deletion of 
an additional particle are no longer equivalent, and sym¬ 
metry of the phase diagram about n/t = 0 and p = ^ 
would be surprising, rather than expected. 

Further qualitative insight into the quantitative be¬ 
haviour revealed in Fig. 2(i) may be obtained by com¬ 
paring the filling fraction as a function oi p/t with that 
for fermions and hard-core bosons, at specific values of 
t±/t. Transects performed at tj_/t = 1 and t±/t = 3, 
as in Ref. 9, are enlightening and are shown in Fig. 2(ii) 
and (iii). 

Consider first Fig. 2 (iii), t±/t = 3, where the behaviour 
of Z 3 anyons at filling fractions 0 < p < ^ is seen to 
closely resemble that of fermions, though with a slight 
deviation towards higher filling fractions in the first con¬ 
ducting phase. As with fermions the repeated exchange 
of Z 3 anyons leads to destructive interference, although 
a sum over three terms rather than two is required in or¬ 
der for the contributing diagrams to sum to zero. During 
real-time evolution of a hard-core gas in the dilute limit, 
particles are relatively free to braid around each other, 
and two particles in sufficiently close proximity to un¬ 
dergo pair exchange may expect their future light cones 
to overlap for a significant interval before encountering a 
third particle. The longer this interval, the greater the 
relative extent of the overlap: 



For a sufficiently dilute gas the time to encounter an¬ 
other particle approaches infinity, under which circum¬ 
stances the overlap of the future light cones approaches 
unity, and the number of times modulo 3 that a given 
pair of anyons undergoes exchange is equally likely to 
be 0, 1, or 2. Full destructive interference, leading to 
precise equivalence of filling fractions for fermions and 
Z 3 anyons, is therefore approached in the limit of a suf- 
hciently dilute anyon gas [57]. At finite but small filling 
fractions, larger numbers of exchanges become less likely 
than smaller numbers for any given pair encounter, so 
P( 2 ) < P(l), and a small residual term arises from zero- 
and single exchanges not matched by double exchanges. 
It is anticipated that for increasing anyon species q in the 
set of models {Z^} and fixed filling fraction p this residual 
term will increase, with the limit g —> 00 corresponding 
to hard-core bosons where pair exchanges are not asso¬ 
ciated with any destructive interference. Thus the Zg 
anyons are anticipated to interpolate between fermions 
and hard-core bosons in the left half of Fig. 2 (iii). 

As with both bosons and fermions, a plateau is ob¬ 
served at p = 5 . Although the behaviour of Z 3 anyons 
still strongly resembles that of fermions in this regime, 
this plateau must be considered a rung-Mott insulator as 


the addition of further particles is not inhibited by lack 
of available states in the band structure, but rather by 
the combination of infinite on-site repulsive interaction 
and the energy penalty associated with obstruction of 
the hopping terms in the Hamiltonian. At the \ow-p/t 
end of the plateau, the anyons’ pseudo-fermionic charac¬ 
ter predominates. However, at the high-/r/f end of the 
p = 1/2 plateau the insulating phase is seen to be more 
stable than its bosonic and fermionic counterparts, as is 
the plateau at p = 1 . 

In this higher filling regime, ^ < p < 1, several effects 
apply. First, the introduction of additional particles re¬ 
duces the amount of hopping which can take place, and 
the increase in chemical potential must be sufficient to 
overcome the associated energy penalty. However, this 
effect applies equally to all species and so does not ex¬ 
plain the unique behaviour of Z 3 anyons. 

Beyond p = 1/2 it is useful to think of the holes as 
the free excitations on an otherwise full lattice, with a 
rung occupancy of n corresponding to 2 — n holes. For 
fermions particle/hole duality applies, allowing the holes 
to be viewed as the carriers of current above p = 1 / 2 . 
Adding particles is therefore favourable as it decreases 
the amount of destructive interference associated with 
hole exchange. Hard-core bosons also exhibit parti¬ 
cle/hole duality, but removal of holes is neutral with re¬ 
gards to interference so a higher chemical potential is 
required to leave the p = 1/2 plateau. Finally, for Z 3 
anyons the holes may still be thought of as the free ex¬ 
citations at higher p but fully-occupied rungs now carry 
a charge of 2 . In this context even relatively simple hole 
trajectories such as 



can generate destructive interference, opposing transition 
into a hole-mediated conducting regime. Similarly, start¬ 
ing from the p = 1 plateau and decreasing p/t the p = 1 
plateau is likewise seen to be more extensive as destruc¬ 
tive interference again opposes the introduction of holes. 
Noting that Z 2 is a special case for which the saturation 
of a rung with particles yields a total rung charge of 0 , it 
is expected that for all other i.e. Zg| 3 <ggz+, destruc¬ 
tive interference will inhibit departure from the p = 1/2 
and p = 1 plateaus. However, this effect will decrease 
with increasing q as the phase angle which is associ¬ 
ated with braiding decreases, attaining hard-core boson 
behaviour in the limit g —>■ 00 . 

In Fig. 2(ii) the behaviour of Z 3 anyons is again 
fermion-like for p <1/2 with deflection towards the hard¬ 
core boson line due to a reduction in destructive inter¬ 
ference. To the extent that a quasi-fermionic description 
applies, the weakening of destructive interference may be 
viewed as deforming the effective band structures, delay¬ 
ing transition from the one-band superfluid phase (SFl) 
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FIG. 3. (Colour online) Probabilities for different separa¬ 
tions of nearest-neighbour holes in the hole-crystal phase with 
tx/l = 1, at (i) ^,/t = 2.06 and (ii) p/t = 2.28. Labels “opp” 
and “same” indicate whether holes are on opposite or the 
same rungs respectively, (hi) Probabilities for different sep¬ 
arations of nearest-neighbour holes as a function of p/t, in 
superfluid and hole-crystal phases with t±/t = 1. 


to the two-band superfluid phase (SF2) and permitting 
the attainment of higher filling fractions in SFl. 

For t±/t = 1 there is no insulating phase at p = ^ due 
to overlap in the effective band structure, and the be¬ 
haviour of the anyonic model continues to qualitatively 
mimic that of the fermionic model right up to p ~ 0.78. 
The role of destructive interference decreases rapidly, 
however, as a higher particle density provides an obstruc¬ 
tion to repeated braiding, and this causes the Z 3 line to 
shift further away from the fermion line and towards the 
hard-core boson line. 

At approximately p/t = 2.04, p = 0.78 there is a phase 
transition from a conventional superfluid phase with pre¬ 
dominantly zero or one vacancies per rung to a phase 
in which the holes are organised in a regular crystalline 
structure (labelled HC). Discontinuity of the 3-site von 
Neumann entanglement entropy shows that this phase 
transition is first-order [Fig. 2(iv)]. In the HC phase, con¬ 
figurations are favoured where sequential pairs of holes 
are located on opposite legs of the ladder, at separations 
which admit bosonic exchange statistics. In the vicin¬ 
ity of the critical point the typical separation is 2 rungs: 
sss or with these configurations being collectively 
denoted 2 °pp ( 2 ®™® would indicate configurations with 
both holes on the same leg). Further peaks in amplitude 
appear at larger separations as the hole density decreases 
[Fig. 3(i)], consistently favouring separations where hole 
exchange yields bosonic statistics either directly or with 
only a small perturbation, e.g. for 6 °pP: . 

In the vicinity of the critical point, the one-hop per¬ 
turbations about this preferred configuration are denoted 
1 °PP, 3 °PP, and 2®™®. As shown in Fig. 3(ii), the prob¬ 
abilities for these configurations loosely converge in the 
vicinity of the critical point, while two-hop perturbations 
are rarer, and the low probability for two holes to appear 
on the same rung, 0 °pp, indicates that hole exchange is 
relatively uncommon. These results are indicative of a 


phase in which holes are regularly spaced and poorly de¬ 
localised, predominantly due to small fluctuations about 
their preferred separations. Finally, the p = 1 plateau 
extends to lower p/t due to destructive interference. 

Generalising to other anyon species Z^, with increas¬ 
ing q it is anticipated that the SFI/SF2 phase transition 
for p < i will continue to migrate right and up, and oc¬ 
cupancy in the SFl phase will shift further towards the 
hard-core boson line as q increases. The SF2 phase will 
become narrower and the SF1/SF2 and SF2/SF1 phase 
transitions will become coincident in the limit g —>■ 00 , re¬ 
covering the hard-core boson behaviour where only SFl 
behaviour is observed throughout. No lateral shift of 
the SF2/SF1 phase boundary was observed, and thus al¬ 
though the (7 —>■ 00 limit of hard core bosons exhibits 
particle/hole symmetry, it may be that this limit is ap¬ 
proached asymmetrically with the phase boundaries con¬ 
verging at some value p/f > 0. For odd q there is a novel 
hole-crystal phase, in which vacancies form a regularly- 
spaced foam purely as a result of the particle exchange 
statistics. The extent of this phase is expected to de¬ 
crease with increasing q. It is characterised by hole sepa¬ 
rations consistent with bosonic exchange statistics, most 
notably (| -I- 5 )°’’^, implying that the filling fraction at 
which the SFl/HC phase transition takes place must sat¬ 
isfy P > q/{q + ^)- This phase therefore vanishes in the 
limit q ^ 00 . For both even and odd Z^ the width of the 
p = 1 plateau is also expected to decrease with increas¬ 
ing q, as the phase angle associated with braiding 
decreases and destructive interference becomes less po¬ 
tent. It is conjectured that hole-crystal phases like the 
one observed here could potentially also appear in full 2D 
systems, with smallest unit cells for Z3 on the square, tri¬ 
angular, and hexagonal lattices appearing as follows: 



In summary, the phase diagram for hard-core non¬ 
interacting Z3 anyons on the two-leg ladder is qualita¬ 
tively distinct both from that of fermions and of hard¬ 
core bosons [9], with destructive interference giving rise 
to extended insulating phases at filling fractions p = 5 
and p = I. A novel first-order phase transition is identi¬ 
fied at high filling fraction, into a phase characterised by 
a regular arrangement of strongly localised holes. Equiv¬ 
alent phases are expected to appear for all Zql^g^odd. 
Qualitative explanations are proposed for all observed 
properties of the Z3 phase diagram, with extrapolative 
predictions being made for Zq anyon models with q > 3. 
These predictions should be readily amenable to testing 
through further numerical simulations. 

The author thanks Gavin Brennen for beneficial dis¬ 
cussions. The simulations performed in this paper made 
use of the net con package [58]. 
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